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Abstract 
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i Predictions are presented for the oscillation phase in neutrino oscillations following pion de- 

I cay at rest and in flight, muon decay and nuclear /3— decay at rest as well as for muon oscillations 

following pion decay at rest and in flight. The neutrino oscillation phases found disagree with 
^ ' the conventionally used value: (/)i2 = Am L/{2P). The same oscillation phase is found for neu- 

■ trino and muon oscillations following pion decay at rest. Shortcomings of previous treatments 

' of the quantum mechanics of neutrino oscillations are commented on. 
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This letter sketches briefly a derivation of spatially dependent interference effects in 
the detection probability of neutrinos and muons following the weak decays of unstable 
'source' particles: pions, muons and /?— radioactive nuclei. It is based on Feynman's path 
amplitude formulation of quantum mechanics in which the probability of transition 
5^ I from a set of initial states / = Yli h to a set of final states F = Ylm fm is given by the 

relation: 
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where kj, j = l,n are (unobserved) intermediate quantum states. Only the essential 
elements of the calculation and the results are presented here. Full details are given 
elsewhere 0]. 

The application of Eqn(l) to pion decay at rest is illustrated in the ideal experiment 
shown in Fig.l. A 7r+ comes to rest in a stopping target T at time to as recorded by the 
counter Ca (Fig. la)). The pion at rest constitutes the initial state of the path amplitudes. 
In Fig. lb) and Fig.lc) are shown two alternative histories of the stopped pion. In Fig. lb) 
the pion decays at time ti into the neutrino mass eigenstate {ui >, of mass mi, and in 
Fig.lc) into the neutrino mass eigenstate \i>2 >, of mass m2, at the later time t2- If 
mi > m2, then, for a suitable choice of the times ti and t2, interference between the path 
amplitudes corresponding to the different physical processes shown in Fig. lb) and Fig.lc) 
will occur when a neutrino interaction e^p takes place at time to and distance L 
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from the pion decay point (Fig. Id)). The final state of the neutrino interaction event is 
also that of the path amplitudes. For this experiment, the path amplitudes corresponding 
to the two alternative histories of the decaying pion are: 

Ai = < e~p\T\nUe >< UeWi > Dixf - Xi,tD -ti,mi) < Ui\Uf_,> 

< z/^^+|r|7r+ > e-^(*'-*«)D(0,t, -to,^^) i = l,2 (2) 

Here < e^p|T|nz/e > and < z/^/i+|T|7r+ > are the invariant amplitudes of the neutrino 
interaction and pion decay processes respectively, < UeWi > and < z/j|i/^ > are amplitudes 
describing the mixture of the fiavour and mass neutrino eigenstates: 

|z/e > = cos6'|z/i > + sin6'|z/2 > (3) 
|z/^ > = — sin6'|z/i > + cos^^|z/2 > (4) 

(only 2-fiavour mixing is considered), and D is the Lorentz invariant configuration space 
propagator f^, Q of a neutrino or the pion. In the limit of large time-like separations or of 
on-shell particles, appropriate to the experiment shown in Fig.l, D ^ exp[— imAr] where 
m is the pole mass of the particle and At the increment of proper time corresponding 

3 

to the path. In the following the additional functional dependence ~ (m/Ar)2 of D in 
the asymptotic region (leading to solid angle correction factors) is neglected. With these 
approximations: 



D{Ax,At,m) ~ exp[-im^/{Aty - (Ax)^] 
= exp[— imAr] 

= exp[-zA0] (5) 

The phase increments, A0, corresponding to the paths of the neutrinos and the pion in 
the amplitudes Ai are: 

2 2 
rnr 777- 

A0r = m,Ar, = -^At, = -^L (6) 
A0. = m^{ti - to) = m^[tD - to) 

= rn^ito -to) - m^^L + (7) 

where mj, Pi and Vi are the mass, energy, momentum and velocity of the mass 
eigenstate Vi and 

2 2 
ffi — 777, 

Po = — = 29.8MeV (8) 

2m^ ^ ^ 

The neutrinos are assumed to follow classical rectilinear trajectories such that Ati = 
L/vi = EiL/Pi and the time dilatation formula At = ^At = EAr/m has been used in 
Eqn(6). In Eqn(7) the neutrino velocity^ Vi is expressed in terms of the neutrino mass to 
order mf. Using Eqns(5-7) the path amplitudes of Eqn(2) may be written as: 



Ai = < e p\T\nVf, >< UeWi >< ^iWtJ, >< ^fj.fJ'^\T\n^ > 
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"Units with /i/27r = c = 1 are used. 
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where the neutrino times-of-flight t{'' = to — ti have been introduced and 



= m^(L - to + to) 



(10) 



Using now Eqn(l) to calculate the transition probability, and integrating over the detec- 
tion time tfl 0, gives, for the probability to observe the reaction u^n — > e~p at distance 
L from the decay point: 



P{e-p\L) = CN{i^]7r) sin2^cos2 0(l-F^(r^)cos0i2" 



(11) 



where 
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and Eqns(3) and (4) have been used. Here CAr(z/;7r) is an L independent normalisation 
factor and Am^ = m\ — m\. The first and second terms in the second member of Eqn(12) 
are the contributions of A(/)J^ and A^J' to the interference phase. The latter is a factor 
of two larger than in the conventional value ^ 0i2 = Am^L/(2P). A0f gives a 
numerically large {mT^/2Po = 2.34) contribution to the oscillation phase. 

The above calculation is readily repeated for the case of detection of muon decay 
e^VgV^j^ at distance L from the tt"*" decay point. The phase increments analagous 
to Eqns(6) and (7) are: 
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The result found for the time-integrated decay probability is: 



P(e+z/eZ7^|L) = C^(/i;7r)(l -P''(r^)sin2ecos05'f) 



(16) 
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The first and second terms in the second member of Eqn(17) are the contributions of 
A(f)J^^^ and A0f to the interference phase. It is interesting to note that neutrino and 
muon oscillation phases are the same for given values of Am^ and L. For oscillation phases 
(f)i2 = 012'^ = I5 the damping factors of the oscillation term, due to the non- vanishing 
pion lifetime take the values F''{T^) = l-1.58x IQ-^^ and F'^(r^) = l-4.4x 10"^^ This 
damping effect is thus completely negligible in typical neutrino oscillation experiments 
with oscillation phases of order unity. 

Formulae like Eqn(ll) have been derived in a similar manner for neutrino oscil- 
lations z7^ — > z7e following muon decay at rest, detected via the process VeP e~^n and 
Ve — >■ z/g oscillations following nuclear /3— decay, detected via the process VeP e+n. The 
results found, for the time-integrated decay probabilities, are: 



P{e~^n,fi\L) = Cn{i^] fi) sm^ 9 cos^ 9 
P{e+n,f3\L) = Cr,{u-p) 



Am ( nin 

1— ^U7t-^l^ 



sin' » + cos' + 2 siri= » cos^ 9 cos ^ " 1 I ^ 



(19) 
(20) 



In these formulae, is the momentum of the detected Ue, Ep is the total energy release in 
the /3— decay process, and damping corrections due to the finite lifetimes of the decaying 
particles are neglected. 

Finally, in Ref. 0, the cases of neutrino and muon oscillations following the decay in 
flight of ultrarelativistic pions was considered. The oscillation probability formulae are 
the same as Eqns(ll) and (16) respectively, with the phases: 

(/.rr in flight = , , J- (21 

c/.^'f (in flight) = ; , ^ (22) 

[m^ — m'^j'^E^ cos t^^ 

The neutrinos or muons are detected at distance L from the decay point, along the 
direction of flight of the parent pion, with decay angles 9i, and 6*^ respectively. 

Because the parent pion and the daughter muon in the decay process vr — >• /ii/ are un- 
stable particles, their physical masses and Wfj, have distributions depending, through 
Breit-Wigner amplitudes, on their pole masses 171-,^ and and decay widths and F^. 
Energy-momentum conservation in the pion decay process then leads to a distribution of 
path amplitudes with different neutrino momenta, an effect neglected in the above dis- 
cussion. In Ref.[0] corrections resulting from the (coherent) momentum smearing due to 

and (incoherent) smearing due to are calculated using a Gaussian approximation 
for the Breit-Wigner amplitudes. The resulting damping corrections due to to the 
interference terms in Eqns(ll) and (16) are found to be vanishingly small. For Am^ = 
(leV)2 and L = 30m (typical of the LNSD [0] or KARMEN g experiments) the corre- 
sponding damping factors are found to be 1 — 1.3 x 10~^^ (for neutrino oscillations) and 
1 — 3.9 X 10"'^^ (for muon oscillations) 0. The damping corrections from the variation 
of Wn are even smaller. 

Corrections to Eqns(ll) and (16) due to thermal motion of the decaying pion and finite 
target and detector sizes have also been evaluated in Ref.0. For neutrino oscillations 
with Am^ = (leV)^ and L = 30m and a room-temperature target, the damping factor of 
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the interference term is found to be 1 — 6.7 x 10^^*^ and the shift in the oscillation phase 
to be 1.2 X 10~^rad. In summary, all known sources of damping of the interference terms 
in Eqn(lO) and (15) are expected to be completely negligible in any forseeable neutrino 
or muon oscillation experiment. 

The first published calculation of the neutrino oscillation phase gave a result (only 
the contribution of the neutrino propagators was taken into account) in agreement with 
Eqn(6) above. A later calculation [rU| assumed instead that the phase of the neutrino 



propagator evolves as Et, i.e. according to the non-relativistic Schrodinger equation. The 
assumption was also made that the different neutrino mass states have equal momenta. 
This leads to the 'standard' formula for the oscillation phase: 0i2 = Am^L/(2P) ^ 
which has subsequently been used for the analysis of all neutrino oscillation experiments. 
The Lorentz-invariant phase of Eqn(6) ~ m?t/E evidently agrees with the result of 



Ref.[|10l in the non-relativistic limit where E c::^ m, but such a limit is clearly inappropriate 
to describe experiments with ultra-relativistic neutrinos. 

The most important difference in the treatment given in the present paper to previous 
ones that have appeared in the literature is allowing the possibility for the different 
neutrino mass eigenstates to be produced at different times. Only in this way can the 
constraints, of both space-time geometry (the detection event is at a unique space-time 



point), and exact energy- momentum conservation in the decay process [0, be satisfied. 

The other new feature is the inclusion of the important contribution to the oscilla- 
tion phase from the propagator of the decaying particle, a necessary consequence of the 
different production times of the different mass eigenstates. 

Many recent treatments of the quantum mechanics of neutrino oscillations have, fol- 



lowing the suggestion of Ref. |T^, used a 'wave packet' description of the spatial distri- 
butions of the source and detector particles. This gives, for the transition probability of 
the complete production-propagation-detection process: 
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Here, and ipf are 'source' and 'detector' wave packets respectively. The use of such 
an equation leads to the prediction, in conflict with Eqn(l) and the fundamental law 
of superposition in quantum mechanics, that different initial or final state particles will 
interfere with each other 0). As is well known, (for example, in Fermi's 'Golden Rules') 
different initial and final states must be (incoherently) summed over at the level of the 
transition probability, not coherently, at amplitude level, as is done in Eqn(23). In a 
previous covariant analysis of neutrino oscillations 0], the invariant neutrino oscillation 
phase, derived from Eqn(6) above, was found to be multiplied by a factor 1/2 on per- 
forming a convolution, according to to Eqn(23), with a Gaussian space-time wave packet 
for the source particle. In this way the standard result of Ref . [|1^ was recovered. 

It may be remarked that the physical interpretation of 'neutrino oscillations' provided 
by the path amplitude description is different from the conventional one in terms of 
'flavour eigenstates'. In the latter the amplitudes of different flavours in the neutrino are 
supposed to vary harmonically as a function of time. In the amplitudes for the different 
physical processes in the path amplitude treatment there is, instead, no variation of 
lepton flavour in the propagating neutrinos. If the mass eigenstates are represented 
as superpositions of flavour eigenstates using the inverses of Eqns(2) and (3), there is 
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evidently no temporal variation of the lepton flavour composition. Only in the detection 
process itself are the flavour eigenstates projected out, and the interference effect occurs 
that is described as 'neutrino oscillations'. In the case of the observation of the recoil 
muons no such projection occurs, but exactly similar interference effects are predicted 
to occur. As previously emphasised |T^, the 'flavour oscillations' of neutrinos, neutral 



kaons and b-mesons are just special examples of the universal phenomenon of quantum 
mechanical superposition that is the physical basis of Eqn(l). 

The most important practical conclusion of the work presented here is that identical 
information on neutrino masses is give by the observation of either neutrinos or muons 
from pion decay at rest. In view of the possibility of detecting muons simply and with 
essentially 100 % efficiency, in contrast to the tiny observable event rates of neutrino 
interactions, the recent indications for z7^ — > z/g oscillations with Am^ ^(leV)^ following 
/i"*" decay at rest could be easily checked by a search for muon oscillations following 7r+ 
decay at rest. For Am^ ~(leV)^, the first absolute maximum of the interference term in 
Eqn(16) occurs at L ~ 8m. Note that Eqn(16) is valid for any muon detection process. 
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Figure 1: The space-time description of i'^ — >■ oscillations following vr"*" decay at rest, in Feynman's 
formulation of quantum mechanics (see text). 
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